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We obtain sharp criteria for transverse stability and instability of line solitons in the discrete non- 
linear Schrodinger equations on one- and two-dimensional lattices near the anti-continuum limit. On 
a two-dimensional lattice, the fundamental line soliton is proved to be transversely stable (unstable) 
when it bifurcates from the X (F) point of the dispersion surface. On a one-dimensional (stripe) 
lattice, the fundamental line soliton is proved to be transversely unstable for both signs of transverse 
dispersion. If this transverse dispersion has the opposite sign to the discrete dispersion, the insta- 
bility is caused by a resonance between isolated eigenvalues of negative energy and the continuous 
spectrum of positive energy. These results hold for both focusing and defocusing nonlinearities via 
^SJ ] a staggering transformation. When the line soliton is transversely unstable, asymptotic expressions 

for unstable eigenvalues are also derived. These analytical results are compared with numerical 
■ results, and perfect agreement is obtained. 
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I. INTRODUCTION 



O: 

CNj I One-dimensional solitons, when viewed in two spatial dimensions, become line solitons which are uniform along 
the line direction (called transverse direction). Thus an important physical question is the transverse stability of 
line solitons to transverse perturbations. It is well known that in homogeneous media, line solitons in the nonlinear 
p ^ . Schrodinger (NLS) equation and other related wave equations are always transversely unstable [32,] (see also [1, H, Q 
for applications in optics and [13, HI] for reviews) . This instability has been observed in recent optical experiments 
Tl]. Ilil . [l9| . In the presence of a one-dimensional periodic potential, many line solitons are still transversely unstable 
ll. |22| . To suppress this transverse instability, various techniques have been proposed [1,113, Hi, 30, 31). In particular, 
it was shown numerically in [s^, HH that when a one- or two-dimensional periodic potential is included in the 
continuous NLS equation, this transverse instability can be completely eliminated if the line soliton bifurcates from 
certain symmetry points of the dispersion surface. But in the corresponding discrete nonlinear Schrodinger (dNLS) 
equations, discrete line solitons in one-dimensional lattices are numerically found to be still transversely unstable 1311] , 
highlighting the difference between continuous and discrete NLS models. 

In this article, we analytically investigate transverse stability of line solitons in the dNLS equations on two- 

■ dimensional (square) and one-dimensional (stripe) lattices. Near the anti-continuum limit, we derive sharp stability 
-iJ ] criteria for these discrete line solitons. We prove for a two-dimensional lattice that the fundamental line soliton is 

■ transversely stable (unstable) when it bifurcates from the X (F) point of the dispersion surface. For a one-dimensional 
] (stripe) lattice, the fundamental line soliton is proven to be transversely unstable for both signs of transverse dis- 

■ persion. These results hold for both focusing and defocusing dNLS equations via a staggering transformation. For 
ILJ ■ unstable line solitons, their unstable eigenvalues are also derived asymptotically. We also investigate the transverse 

stability of line solitons numerically both near the anti-continuum limit and away from it. Near the anti-continuum 
limit, the numerical results fully agree with the analytical results. Away from the anti-continuum limit, we reveal 
^ ] additional bifurcations of unstable eigenvalues which cannot be captured by the theoretical analysis. 

This article is organized as follows. In Section |lll we consider the transverse stability of discrete line solitons on a 
two-dimensional lattice. We show that the entire solution family bifurcating from the F point is transversely unstable, 
whereas the solution family bifurcating from the X point is transversely stable in the anti-continuum limit. In Section 
mil we consider the transverse stability of discrete line solitons on a one-dimensional lattice, and show that they 
are always unstable. Numerical results and their comparison with the theory are reported in Section IIVI Section IVl 
concludes the paper with discussion of open problems. 

Before we start, we first introduce some mathematical notations which will be used in later analysis. If {V'njnGZ is 
a bi- infinite sequence (i.e., a sequence which is infinite in both directions), and Z is the set of integers, then denotes 
the vector for this sequence in some vector space such as Z^(Z) or, more generally, PCZ) for p > 1. Here denotes 

1/2 

the space of bi-infinite squared-summable sequences with the norm \\tp\\i2 = (X^nez IV'nj^) the inner product 

= X^nez 'Z^"'''"' "^ith the overbar for complex conjugation, and Zp(Z) denotes the space of sequences with the 
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II. TWO-DIMENSIONAL LATTICE 



In this section, we study transverse stabihty of hne sohtons in the dNLS equation on a two-dimensional lattice, 

^ f^(''^77i+l,n ~t- Wm — l,n ^^ ^m,7i+l ^" ^m.n~l ^^m.n) |w7n,n| '^m^n — 0, (1) 

where {m,n) G Z^, Um,n are complex- valued amplitudes that depend on the evolution time i, and e is the lattice- 
coupling constant. Here the sign of nonlinearity has been normalized to be unity through a scaling of e and t. The 
anti-continuum limit e = of zero coupling between lattice sites was found to be very attractive for many analytical 
studies on the existence and stability of discrete solitons in the framework of the dNLS equation 0, O, [25L J2p|. 
Detailed account of mathematical results obtained in the ant i- continuum limit can be found in the monograph [24| . 

In the above dNLS equation, the defocusing case e < can be mapped to the focusing case e > by the staggering 
transformation 

um,n{t) = (-l)''"+"^^™,„(i)e"'"*. (2) 

If u solves the dNLS equation ([T|), then v solves the same equation with e replaced by — e. Thus, in what follows, we 
will consider the focusing case (e > 0) only. 

The linear dispersion surface of the dNLS equation ([1]) is given by the function 



u{k,p) = e(4 - 2cos(fc) - 2cos(p)) = 4e 



sin^ I - 



where wavenumbers {k,p) reside in the first Brillouin zone [— 7r,7r] x [— tt, tt]. This dispersion relation can be derived 
by substituting the discrete Fourier modes Um,n{t) — e*'""+*P"^*'^* into the linear dNLS equation ([T|). 

To understand bifurcations of stationary line solitons in Eq. ([T]), we need to classify the stationary points of 
the dispersion surface, where yuj{k,p) = 0. In the semi-open Brillouin zone (— 7r,7r] x (— 7r,7r], there are only four 
stationary points, which are commonly labeled as F, X, X' , and M . 

F: {k,p) = (0, 0) is the minimum point of the dispersion surface with ^(0, 0) — 0; 

X: {k,p) — (0, tt) is a saddle point of the dispersion surface with w(0, tt) = 4e; 

X': {k,p) — (tt, 0) is the other saddle point of the dispersion surface with a;(7r, 0) — 4e; 

M: {k,p) — (7r,7r) is the maximum point of the dispersion surface with w(7r,7r) = 8e. 

Discrete line solitons may bifurcate from any stationary point provided that the effective continuous NLS equation 
is focusing Section 1.1.2]. Let us consider each of the possibilities. For definiteness, we assume that the line 
soliton is localized along the m-direction and uniform along the n-direction. 

F: For (fc,p) = (0,0), we substitute = e'^ *V'm and obtain the stationary ID dNLS equation 

- fi'^tpm + eilpm+l + Tpm-l - 2V'm) + \^m\'^1pm = 0, (3) 

which admits discrete solitons for any e > and ^ 7^ [13, IS^- Moreover, for fixed e > 0, the fundamental 
discrete soliton is approximated by the NLS soliton 



'i/'m V2/i sech [ ) as fi 0. (4) 



This approximation was rigorously justified in the recent work [4| (see also [23, Section 2.3.2]). 

X: For (fc,p) — (0,7r), we substitute Um.n{t) = (— l)"e*'^'^ ~^'''*'!/'ni and obtain the same stationary dNLS equation 
([3]), which admits the discrete solitons for any e > and 7^ 0. 

X': For {k,p) = (tt, 0), we substitute u„ „(<) = (— l)™e*(^^^^^^^*-0m and obtain the stationary ID dNLS equation 

H'^lpm - + tpm-l - 2^m) + \tpm\'^'4'm. = 0. (5) 

This stationary equation admits no discrete solitons for any e > [23, Lemma 3.10]. Indeed, by projecting Eq. 
^ to Tp and denoting {A^)m = ipm+i + ''Pm-i — 2'(/'m, we obtain a contradiction 

/i2||^||?.-t-e(V,(-A)V) + ||V||f4=0, 

where each term on the left side is positive definite. 
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M: For {k,p) = (tTjTt), we substitute Um,n{t) = (— l)'"+"e'( ^'^•'*V'm and obtain the same stationary dNLS 
equation ([5]), which admits no discrete sohtons for any e > 0. 

From the above analysis, we see that only two bifurcations of fundamental discrete line solitons occur and the 
bifurcation points are F and X. In the absence of transverse perturbations, these fundamental line solitons are stable. 
In the following, we will analyze transverse stability of these fundamental line solitons in the ant i- continuum limit 
e — >■ for fixed /j, > (or equivalently, /i — >■ oo for fixed e > 0). 

Before the transverse-stability analysis in the ant i- continuum limit, it is useful to recall the transverse-stability 
results in the opposite (continuum) limit that arises when e — > oo for fixed /i > (or equivalently, /i — > for fixed 
e>0). 



F: For (fc,p) = (0,0), we substitute 



into Eq. ([T]). Assuming smoothness of the envelope function U{X,Y,t)^ we obtain an elliptic 2D NLS equation 

for U{X, Y, t) a.s e oo: 

.dU d'^U d'^U 2-.rr r. 

The line soliton (U) is transversely unstable in this elliptic NLS equation © due to neck-type instability (see 
[itI [29} and references therein) . 

X: For {k,p) = (0, tt), we substitute 

umAt) = (-i)"c/(x,y,r)e^(^'-4.)t^ x^^,Y = ^ 



into Eq. ([T]). Assuming smoothness of the envelope function U{X, Y, t), we obtain a hyperbolic 2D NLS equation 
for U{X, Y, t) SiS e ^ 00: 

.dU d'^U d^U n 

The line soliton dH) is also transversely unstable in this hyperbolic NLS equation (O due to snaking-type 
instability (see [lO, 29] and references therein). 

From reductions to 2D NLS equations © and ([7]), we see that discrete line solitons (|3]) are always transversely 
unstable in the continuum limit. Thus it is surprising that discrete line solitons were reported to be transversely 
stable far from the continuum limit when they bifurcate from the X point of the dispersion surface [soj . Line solitons 
bifurcated from the F point, however, remain transversely unstable for all values of e (i.e., both near the continuum 
limit and away from it) [30] . Below we shall prove these numerical observations by rigorous spectral-stability analysis 
that relies on the count of eigenvalues of negative energy @, [l^ [2^ . In addition, asymptotic expressions for unstable 
eigenvalues will also be derived in the anti-continuum limit. 



A. Instability of line solitons bifurcating from the F point 

Discrete line solitons bifurcating from the F point are of the form 

Um,„(t) =e*'^'Vm, (8) 

where satisfies the stationary ID dNLS equation (jS]). It can be easily shown that {V'mlmez in these discrete solitons 
is real-valued (up to multiplication by e*" for real a, i.e., a € K) [13, Lemma 3.11]. Perturbing these line solitons as 

Uni.nit) = e''"'* [ipm + 1'm,«(t)] , 

and substituting it into the dNLS equation ([T}, we obtain the linearized dNLS equation as 

.dUm^ - IJ.'^Vm,n + e(Wm+l,n + Wm-l.n + Vm,n+l + Vm,n-l - ^Vm,n) + V'm(2^'m,n + Vm,n) = 0. 
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For normal modes 

we obtain the standard form of the eigenvalue problem 

L+ip)U = -XW, L^{p)W ^ XU, (10) 
where L±{p) are p-dependent ID discrete Schrodinger operators, 

{L+{p)U)m = -e [Um+l + Um-1 + 2 COs{p)Um - 4f7„] + fM^Um - i^ljJm, , x 

{L_{p)W)m = -e [W^+i + Wra-i + 2 cos(p)W^™ - AWm] + iJi'^W-^ - ^IW^. ^^^> 

It is easy to see that eigenvalues A in the above linear-stability problem always appear as quadruples (A, A, — A, — A) 
when A is complex or as pairs (A, —A) when A is real or purely imaginary. 

Among the two parameters /i and e in the above eigenvalue problem, the ratio e//x^ is invariant with respect to a 
scaling transformation. The anti-continuum limit corresponds to the limit of e//i^ — 0. Without loss of generality, 
we fix /i = 1 and consider small values of e > below. 

We are interested in transverse stability of the fundamental line soliton ■(/;,„, which is positive for all m € Z and 
confined to a single lattice site, say at m = 0, in the anti-continuum limit e 0. Because the stationary equation 
([3]) is analytic in e and polynomial in ip, whereas the difference operator is bounded, the dependence of '0 on e is real 
analytic near e = (2^ . Theorem 3.8]. Using the regular perturbation method, we can easily obtain the power series 
expansion for ip as 

1pm = Srn,a + ^('^m,! + (^m.O + Sm.-l) + 0{e'^), (12) 

where 5m,m' is the Kronecker notation with Sm,m' = 1 for m = m' and otherwise. 

We shall now present the instability theorem for fundamental discrete line solitons bifurcating from the F point. 
This fundamental line soliton exists for any e > [3,[23| (see also [3 Theorem 3.12]). Our instability theorem below 
applies to all values of e > 0, except that the asymptotic expression for the unstable eigenvalue is valid only near the 
anti-continuum limit e — )• 0. 

Theorem 1. Consider the fundamental discrete line soliton f^i bifurcating from the F point in the dNLS equation 
For any e > 0, there is po{e) € (0,7r] such that for any p e (— po(e),Po(e))\{0} the linear- stability problem mO\) admits 
a symmetric pair of real eigenvalues ±A(e,p) with A(e,p) > 0. Hence this fundamental line soliton is transversely 
unstable for all e > 0. In addition, po{e) = tt z/0 < e < i. Furthermore, for any p G [— 7r,7r], the eigenvalue X{e,p) 
has the following asymptotic expansion in the anti- continuum limit, 

X^{e,p) = Sesin^ (|) -f 0{e^) as e ^ 0. (13) 

Proof. We first rewrite operators L±{p) in (fTTj) as 

i±(p) -i±(0) + 2e[l-cos(p)]. 

These are bounded operators from Z^(Z) to P(Z), which have both continuous and discrete spectra. 

The stationary equation ([3]) is simply L_(0)V' = 0. Because ip is positive, is at the bottom of spectrum of i_(0), 
so that i_(0) is non-negative [13, El]- By the perturbation theory, L^{p) is strictly positive for any p e [— tt, 7r]\{0} 
and e > 0. On the other hand, i+(0) has at least one negative eigenvalue because 

(L+(0)V^,^) =-2||V||f4 <0, 

where ||'0|lf4 — X^nez IV'nl"'- Moreover, in the limit e — 0, only one negative eigenvalue of i+(0) exists, which is the 
eigenvalue —2 associated with the central site m — Q. By the variational arguments [is'l, this negative eigenvalue 
persists and remains the only negative eigenvalue of i+(0) for any e > 0. Since L^{p) > L+(0), L-^-{p) has at most 
one negative eigenvalue and no zero eigenvalues. It follows from the stationary equation ^ with fJ- = 1 that 

where A is the ID discrete Laplacian. Thus we obtain 

(i+(p)V, V-) = -2||Vllf4 + 2e [1 - cos(p)] < 2 {e [1 - cos(p)] - 1} UWf,, 
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hence L+{p) admits a negative eigenvalue for any p €E [— 7r,7r] if < e < ^ and for at least small p if e > is 
arbitrary. In other words, for any e > 0, there is Po{() G (0,7r] such that L^{p) has a negative eigenvalue for any 
p e (— po(e)7Po(e))- Moreover, po(e) = tt at least for < e < i. 

For p — Q, the linear eigenvalue problem (jlOp admits zero eigenvalue of algebraic multiplicity two for any e > 
because L_(0)V' = and 

This zero eigenvalue is destroyed when p and this may cause instability when splitting of this double zero 
eigenvalue occurs along the real axis. Using the negative index theory [1, [3; 0, 113] j we obtain: 

where N^^^ (-^reai) ^^'"^ ^^"^ numbers of real positive eigenvalues A with positive (negative) quadratic form {L^{p)U, U) 
at the eigenvector (U,W) of the eigenvalue problem ()10p . N-^^^ is the number of purely imaginary eigenvalues A 
with Im(A) > and negative quadratic form (L+(p)[/, [/), and iVcomp is the number of complex eigenvalues A with 
Re(A) > and Im(A) > 0, counting their algebraic multiplicities. Note that eigenvalues contributing to N^^^ are 
called eigenvalues with a negative Krein signature [l6| . The eigenvalue-counting formula p4p follows directly from 
Theorem 4.5 of [24i] because operators L±{p) have no zero eigenvalues for any p 7^ 0. 

The preceding computations show that n{L_{p)) — 0, and n{L^{p)) = 1 for p e Po(e)jPo(£))- In these cases, the 
index formula (ITi|) yields 

^.;al = 1> ^rtal = ^imag = ^^comp = 0, 

which proves the statement of Theorem 1 on transverse instability. It remains to justify the asymptotic expansion 
([T3| for the real positive eigenvalue A(e,p) as e — > 0. 

When e = 0, the spectral problem (|10p with p, = 1 has three points in the spectrum: A = of algebraic multiplicity 
two and A = ±i of infinite algebraic multiplicity. Continuous spectral bands bifurcate from the points A = ±i for e 7^ 0. 
This bifurcation was studied in detail in the recent work [26], and no unstable eigenvalues arise in this bifurcation. 
We shall now calculate the splitting of the double zero eigenvalue for any fixed p > and small e > 0, using the 
expansion (|12l) near the ant i- continuum limit. 

We rewrite the eigenvalue problem (fTO|) with /i = 1 at the central site m = as follows: 



2C/o + e[Ui + U-i + 2 cos(p)f/o + 2Uo] + 0{e^)Uo = XWq, 
e [Wi + W-i + 2 cos(p)Wo - 2Wo] + C(e^)VFo ^ -XUq. 

By using the scaling transformation U = ^/eU, W = W, and A = -\/e A, the system can be rewritten in the equivalent 
form: 

2Uo + e{Ui + W_i + 2 cos{p)Uq + 2Wo) + 0{€'^)Uq = AWq, 
Wi + W-i + 2cos(p)Wo - 2Wo + C'(e)Wo - -MAq. 

At the adjacent sites m = ±1, the linear eigenvalue problem pO)) is 

U±i - e [U±2 +U0 + 2 cos{p)U±i ~ AU±i] + 0{e^)U±i = -AW±i, 
W±i - e [W±2 + Wo + 2 cos(p)W±i - 4W±i] + 0{e^)W±i = eAU±i, 

since = Similar equations can be written for any m ^ 0. 

For A = 0(1), we have the reduction ll±m = 0{e"^)Uo and >V±m = 0(e™)Wo for any m G N, which enables us to 
close the leading-order equations for {Uq, Wo): 

2Uo + 0{e)Uo = AWo, 

[2cos(p) - 2] Wo + 0(6) Wo = -AZ^o- 

After eliminating Wo, we obtain the algebraic equation for A as 

A^ = 2(2 - 2 cos(p)) + 0(e) = 8 sin^ (|) + 0(e), 
which then yields the asymptotic expansion (1131) . □ 
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Remark 1. For any fixed p € [— tt, 7r]\{0}, we obtain > as e — ^ 0, which guarantees spectral instability of these 
discrete line solitons for small e > 0. Note that the asymptotic formula is not uniform as e — >■ and p — >■ 0, and 
a different perturbation theory is needed in the limit p — )■ for fixed e > (see and references therein). 

Remark 2. Using the resolvent analysis from [26. 1. one can show that the continuous spectral bands are the two 
line-segments on the imaginary axis, 

iX e [-1 - 2e(3 - cos(p)), -1 - 2e(l - cos(p))] U [1 + 2e(l - cos(p)), 1 + 2e(3 - cos(p))], 

whereas no discrete (isolated) eigenvalues bifurcate out from the point X — ±i as c ^ 0. 

Remark 3. In the continuous limit e — > +oo, the discrete line solitons are asymptotically described by the elliptic NLS 
equation [SD, where unstable real eigenvalues are restricted to the p-interval (— po(e),Po(e)\{0} with po{e — > oo) = \/3 
(for fJ. = 1) fl^ . [Hy. Therefore, po{e) < tt for sufficiently large positive e. 



B. Stability of line solitons bifurcating from the X point 

Discrete line solitons bifurcating from the X point arc of the form 

U„,.n{t) = (-l)"e'('^'-4')*Vr„, (15) 

where is a, real-valued solution of the stationary ID dNLS equation ([3])- Notice that these solitons at adjacent lattice 
sites along the transverse n-direction are out-of-phase with each other, which contrasts the line solitons bifurcating 
from the T point, where the solitons at adjacent lattice sites along the transverse direction are in-phase with each 
other. 

Linearizing the dNLS equation ([ij around this solution, we substitute 
and obtain the linearized dNLS equation 

For normal modes ([9]), we obtain the eigenvalue problem 

L+{p)U^~XW, L_{p)W^XU, (16) 

where 

{L+{p)U)m = -e [Um+l + Um-l - 2cOs(p)C/„] + Ai^C/m - SV-^t^m, 

{L^(p)W)m = ~e [Wm+l + " 2cOs(p)VK„] + ^i^'Wm " VmWm- ^ ' 

Again, we are interested in transverse stability of the fundamental line soliton tp which is positive and given by the 
power series expansion (1121) in the anti-continuum limit e — >■ 0. As before, we will fix = 1 without loss of generality. 
The next theorem guarantees stability of this fundamental discrete line soliton for small values of e. 

Theorem 2. Consider the fundamental discrete line soliton kl5) bifurcating from the X point in the dNLS equation 
(Qp. There exists eo > such that for any e € (0,eo) and p £ [tt, tt], the linear- stability problem JJQj) does not admit 
any unstable eigenvalues, thus the fundamental line soliton for small values of e is transversely stable. This stable line 
soliton possesses a pair of discrete imaginary eigenvalues ±za;(e,p) of negative Krein signature. Moreover, for any 
p G [— 7r,7r] and small e, this eigenvalue u)[e,p) has the following asymptotic expression. 



w^(e,]5) = 8esin^ (^IJ +0(e^) as e ^ 0. (18) 

Proof. We first rewrite operators L±{p) in p7)) as 

L±{p) = L±{Q)~2e [l-cos(p)]. 

Because I/_ (0)1/1 = and ip is positive, is the lowest eigenvalue of i-(0) for any e > [l3,[lB|- By perturbation 
theory, i-(p) has exactly one negative eigenvalue for any p G [— 7r,7r]\{0} and small positive e. On the other hand. 
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since and i+(0) for the line soliton (|T5|) are the same as those for the hne sohton ([8|), the variational arguments 
from 13] imply that i+(0) has exactly one negative eigenvalue and no zero eigenvalue for any e > 0. Therefore, 
L+(p) has a single negative eigenvalue for any p e [— 7r,7r] and small positive e. 

For any p ^ 0, we again use the eigenvalue-counting formula (jl4p . which equally applies to the linear eigenvalue 
problem ([16]). The preceding computation shows that there is eo > such that n{L^{p)) — 1 and n{L^{p)) ~ 1 for 
any e G (0, eo) and p G [— tt, 7r]\{0}. Since eigenvalues in the spectral problem (fT6|) appear as quadruples (A, A, —A, —A) 
for complex A^ and as pairs ±A for real A^ and since the zero eigenvalue for p = has algebraic multiplicity two, this 
zero eigenvalue splits along the real or imaginary axis as a pair of simple eigenvalues for p ^ 0. Combining this with 
the eigenvalue-counting formula (I14p . we easily see that this splitting occurs along the imaginary axis, and for any 
ee (0,eo) andpe [-tt, 7r]\{0}, 

which proves the transverse-stability statement in Theorem 2. Note that the imaginary eigenvalues of negative Krein 
signature persist on the imaginary axis, unless they coalesce with other eigenvalues of positive Krein signature or 
continuous spectral bands. 

Next we prove the asymptotic expansion ()18p for the imaginary eigenvalue iuj{e,p) as e — ^ 0. When e = 0, the 
spectral problem (llOp with n = 1 has three points in the spectrum: A = of algebraic multiplicity two and A = ±i 
of infinite algebraic multiplicity. For small e, we only need to compute the splitting of the double zero eigenvalue for 
any fixed p € [— 7r,7r], using the expansion (jl2p near the anti-continuum limit. 

Repeating the perturbation expansions and using the scaling transformation U = y/eU, W — W, and A = -^/e A, 
we obtain the linear eigenvalue problem at the central site to = 0: 

2Uo + e[Ui+U^i-2 cos{p)Uo + O^o] + 0{e^)Uo = AWq, 
Wi + W-i - 2 cos(p)Wo + 2Wo + 0{e)Wo = -AUq. 

Similar to the proof of TheoremlU for A = 0(1), we have the reduction U±m = 0{e"'-)Uo and W±m = 0{e"^)Wo for 
any m G N, hence the above equations yield 

A^ = -2 [2 - 2 cos(p)] + 0{e) ^ -8 sin^ (|) + 0(e), 

which yields the asymptotic expansion ()18|) . □ 

Remark 4. For any p G [— 7r,7r]\{0} and small values of e, we get A^ < 0, which gives imaginary eigenvalues 
±iuj{e,p). It is easy to see that these imaginary eigenvalues have negative Krein signature, meaning that the quadratic 
form {L^{p)U, U) at the eigenvector ([/, W) is negative. For small values of e, these imaginary eigenvalues are hounded 
away from the continuous spectrum bifurcating out of the points ±i, which guarantees spectral stability of discrete line 
solitons for small positive e. 

Remark 5. Using the resolvent analysis from fEdl l. one can show that the continuous spectral bands are located at the 
two segments on the imaginary axis: 

iX G [-1 - 2e(l + cos(p)), -1 + 2e(l - cos(p))] U [1 - 2e(l - cos(j3)), 1 + 2e(l + cos(p))], 

and no discrete (isolated) eigenvalues bifurcate out from the points A = ±i as e 0. 

Remark 6. In the continuum limit e — +oo, discrete line solitons U5]) from the X point in Eq. (Q)) are asymp- 
totically described by the hyperbolic NLS equation where line solitons are transversely unstable for any nonzero 
transverse wave number p fid . l2^J . Hence unstable eigenvalues must appear for these discrete line solitons at suffi- 
ciently large positive e. These unstable eigenvalues can appear through collisions of imaginary eigenvalues of negative 
Krein signature with the continuous spectral band or with additional imaginary eigenvalues of positive Krein signature. 



III. ONE-DIMENSIONAL (STRIPE) LATTICE 

In this section, we consider transverse stability of line solitons in the dNLS equation on a one-dimensional lattice 
with continuous transverse dispersion. The mathematical model for this problem is 

+ e(Mm+l + Ura-1 - 2Um) + K + = 0; (19) 
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where m € Z, the complex variable Um depends on the evolution time t and the transverse coordinate y. Here the 
sign of nonlinearity has been normalized to be unity through a scaling of e, k and t. By the staggering transformation 

u,n{y,t) = i-iy''v^iy,t)e~^''\ (20) 

we can map the dNLS equation (|19|) for u with e < to the same equation for v with e > 0. Thus we set e > below 
but consider both positive and negative values of the transverse dispersion parameter k. Through a scaling of y, we 
normalize k so that k — ±1. 

Transverse instability of line solitons was reported in "sT] for k = — 1 and any e > 0. We shall prove this numerical 
observation by rigorous study of spectral stability. We shall also study the case k ~ +1 for completeness. 

First, by inserting the discrete Fourier modes Ura{y,t) = e*'=™-*'^* into the linear dNLS equation p^ . we find that 
the discrete-dispersion relation is 

= 2e [l-cos(fc)], 

where the wavenumber k is in the first Brillouin zone k e [— tt, tt]. For e > 0, discrete line solitons bifurcate from the 
minimum of this dispersion curve towards negative values of uj. Therefore the discrete line solitons are of the form 

Um(y,t) =e''^'Vm, (21) 

where ^ is a real- valued solution of the stationary ID dNLS equation ([3]). Linearizing around this solution, we 
substitute 

UrniV, t) = e''''* [-0m + Vm{y, t)] 

into the dNLS equation ([T9|) and obtain the linearized dNLS equation 



■ dVm 2 I / I o ^ I ^^^m , / 2 /o , - \ n 
fJ- Vm+ e{Vm+l + Vm-1 " 2u™) + K + tp^{2Vrn + Wm) = 0. 

ot oy^ 

For the normal mode 

we obtain the linear-stability eigenvalue problem 

L+ip)U = -XW, L^{p)W = XU, (22) 



where 



{L+{p)U),n = -e(C/m+l + C/m-1 - ^Um) + (M^ + KP^Wm " 30^C/™, 

{L-{p)WU = -e(W^m+i + Wrn-1 - 2Wm) + (a*' + KJ5')VK™ - ^^l^Wm- 



(23) 



As before, we set /i = 1 by variable rescaling and consider the fundamental line soliton represented by the power series 
expansion (fT2|) for small e. 

When e = 0, the eigenvalue problem (|22p with /i = 1 has four points in the spectrum: two simple eigenvalues at 
A = -izyj Kp^{2 — np"^) and two other eigenvalues of infinite algebraic multiplicities at A = ±i(l -I- Kp^). 

If K = 1, the simple eigenvalues A = ±p-y/2 — p^ are real for < < 2, thus the discrete line soliton ([2T|) is 
transversely unstable even in the unperturbed (e = 0) case. These real eigenvalues persist for small e. The following 
theorem shows that the transverse instability of discrete line solitons pT|) with k = 1 holds for any e > 0. 



Theorem 3. Consider the fundamental discrete line soliton \21\) in the dNLS equation H9\) with k = 1. For any 
e > 0, there is po{e) > such that for any p S {—po{e) , Po{e))\{0} the linear- stability problem i2'^) admits a pair of 
real eigenvalues ±A(e,p) with X{e,p) > 0, thus this line soliton is transversely unstable. In addition, for small e, po{^) 
and A(e,p) are given asymptotically by 



Po{e) ^V2 + 0{e), X{e,p) - p^/2 - p^ + 0{e) as e ^ 0. (24) 
Proof. We first rewrite operators L±{p) in (^5)) as 



L±{p)^L±{0)+p\ 
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Similar to the proof of Theorem [1] we can see that L-{p) is strictly positive for any p and e > 0. On the 
other hand, for any e > 0, i+(0) has a single negative eigenvalue — /?, where (3 > 0. In particular, when e = 0, the 
negative eigenvalue with /? = 2 is associated with the central site m — 0. Thus by denoting po(e) = see that 

L+{p) has exactly one negative eigenvalue for p £ (— po(e),Po(e)) and is strictly positive for any \p\ > po(e)- The 
eigenvalue-counting formula ([T4)) then yields that when p S (— |5o(e),Po(e))\{0}, 

^.^al - h ^+al = ^imag = ^comp = 0. 

The asymptotic expansions p4p directly follow from the preceding computations at e = and the analyticity of the 
linear eigenvalue problem (j22p in e. □ 



Remark 7. The asymptotic expansion X{e,p) = p\Jl — p^ + 0[e) works equally well for \p\ > po{e), where this 
X{e,p) is purely imaginary. These imaginary eigenvalues have positive Krein signature and are hounded away from 
the continuous spectrum located at 

i\e [-(l+p2 + 4e),-(l+p2)]u [l+p^l+p2+4g]. 



If K = —1 and e = 0, the simple eigenvalues A — ±ipy^2 + jp are purely imaginary, so are the eigenvalues 
A — ±i(l — p^) of infinite algebraic multiplicities. These eigenvalue branches intersect at p = ±Pc, where Pc = |. 
When p ^ ±pc and < e ^ 1, the simple eigenvalues persist on «R, whereas two continuous spectral bands bifurcate 
from the non-simple eigenvalues A = ±i(l — p^) along the two segments on as 

iXe [l-p^l-p2+4e]U [-(l-p2+4e),-(l-p2)]^ 

However, when p = ±pc and < e 1, a resonance occurs between these simple and non-simple eigenvalues, and 
as a consequence, complex (unstable) eigenvalues bifurcate out. Notice that the simple eigenvalues A = ±ipy/2 + p'^ 
have negative Krein signature, whereas the non-simple eigenvalues A = ±?(1 — p'^) have positive Krein signature. 
This bifurcation of complex eigenvalues due to collision of eigenvalues with opposite Krein signatures is a common 
phenomenon in Hamiltonian systems [l^, [2^ . 

The following theorem guarantees instability of discrete line solitons ((2T|) in the dNLS equation (|T9)) with k = — 1 
for small values of e > 0. This instability is caused by complex eigenvalues with small real parts, and it occurs for 
intermediate values of transverse wavenumbers p. 

Theorem 4. Consider the fundamental discrete line soliton \21\) in the dNLS equation \19\) with k = — 1. There 
exists Co > such that for any e G (0, eo), there exist p^{e) with ordering < p^(e) < pt i^) < +oo, so that for any 
\p\ G (p~(e),p^(e)) the linear- stability problem admits a quartet of complex eigenvalues ±A(e,p), ±A(e,p) with 

RcA(e,p) > and ImA(e,p) > 0. In addition, when e — > 0, p^{e) and \{e,p) are given asymptotically by 




and 



p^{e) = ^ + -\l±^\+0{e% (25) 



2e 



A(e,p) = -i + -(14 + 17^) + -v/15- 4(1-2(5)2 + 0(^2), (26) 
4 15 15 

where 6 = e^^{p^ — j) = 0(1). Furthermore, the most unstable eigenvalue Xmaxi^) occurs at the transverse wavenum- 
bers ±pmax{(-), where Xmax{e) andpmax{e) are given by 

Xmaxie) = J* + e (-^ + ^i] + 0{e^), p,nax = ^ + + 0{e^). (27) 



Proof. Modifying the arguments from the proof of Theorem [3J we have now 

L±ip) ^ L±iO) - p\ 

Therefore, for sufhciently small e, there is po{e) > such that the operators L±(jj) have exactly one negative eigenvalue 
for ah p e (-po(e),Po(e))\{0}. Note that po{e) = I -h 0{e) as e ^ 0. 
The eigenvalue-counting formula (|14p yields now 

^^imag + ^comp = 1, N+^, = TV^j = 0, p G {^po{e) , p,{e)\{0} . 
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The preceding computations and the analyticity of the linear eigenvalue problem (p2)) in e imply that for sufficiently 
small e, there are pf{e) with ordering < p~(e) < _p^(e) < Po{^) such that 

^i7nag = 1' ^comp = 0, for \p\ € (0,p-(e)) and (p+(e),po(e)), 

where p^{e) = ^ -\- 0{f) as e 0. For these values of e and p, the discrete line solitons are spectrally stable. It 
remains to show that 

^mag = and A^comp = 1 for IpI e (p~(e),p+(e)) 

due to a resonance between eigenvalues of negative and positive Krein signatures. 
First we introduce a scaling transformation 

P = - + ed, \ = — + lej, U,n = , W,n = 



where (5, 7 = 0(1). Under this transformation, the eigenvalue problem (|22|) for fJ. = I and k = — 1 becomes 

- e{am+i - 2am + flm-i) - (1 + 2e)(5„i,o(2ao + 60) + 0{e^){2a„i + 6m) = e(7 + <5)am, (28) 
-e(6m+i - 2bra + 6m-i) - (1 + 2e)5m,o(ao + 26o) + C(e^)(am + 26m) = - + £7 - £'5^ ^m- (29) 

From the second equation (1291) . we obtain 

60 = -2(1 - 2e + 2e7 - 2e5 + O{e'^))ao, (30) 

whereas 6±m = C'(e™)6o for any m e N. The first equation (|28p for any m 7^ produces the second-order difference 
equation 

-(am+i - 2am + flm-i) + 0{e)am = (7 + <^)am, m G Z\{0}, 
which admits a unique decaying solution for both to — > 00 and to — > — 00: 

am ^aoe-"!"!, to e Z\{0}, 
where p is a unique root of the transcendent equation 

-f + 5 = 2-2 cosh(p), Rc{p) > 0. (31) 
To obtain the value for p we close the first equation (1^51) at to = 0: 

-26(6-" - l)ao - (1 + 2e + 0{e^)){2aQ + 60) = e(7 + <5)ao. 
Utihzing ([30)) . this equation becomes 

-2{e-P - 1) - 4(1 - 7 + 5) + 0(e) = 7 + 5. 
Substituting (1311) and neglecting the 0(6) term, we convert this equation to a quadratic equation for z = e^: 

3z2 + 4(25- l)z + 5 = 0, 

which admits two possible solutions 



2(1 - 26) ± - 4(1 - 2(5)2 

3 ■ 

With the help of (I5T]) . these solutions produce expressions for 7 as 



14 + 17(5 + z\/l5- 4(1 - 2(5)2 
^ = ^ ' 

which are complex- valued if (1 — 2(5)2 ^ These expressions yield the asymptotic approximations ([25]) and (|26|). 
From (|26p . we see that the most unstable eigenvalue occurs at (5 = ^, which yields the asymptotic approximation 
(EZl). □ 
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IV. NUMERICAL RESULTS 



In this section, we present numerical results on transverse-stability eigenvalues of discrete line solitons in one- and 
two-dimensional lattices for various values of lattice coupling parameter e (with fixed /i = 1). These numerical results 
are shown to be in good agreement with the analytical results both qualitatively and quantitatively. 



A. Numerical results for the dNLS equation on a two-dimensional lattice 



First we consider discrete line solitons ([5]) bifurcating from the F point in the dNLS equation ([T]). At three values 
of e, eigenvalues of the spectral stability problem (jlOp for various transverse wavenumbers p in the interval [0,7r] are 
presented in Fig. [T] (eigenvalues for negative p are the same as those for positive p). We see that when e — 0.1, a 
single pair of real eigenvalues exist for all values of p in (0,7r], in agreement with Theorem [1] These real eigenvalues 
closely match the asymptotic formula (IT51) in Theorem [T] (middle left panel). When e = 1, this pair of real eigenvalues 
exist only in the interval of < p < po: where po ~ 2.51. For p > po, these real eigenvalues become purely imaginary. 
When e = 4, the p-interval of real eigenvalues further shrinks to (0,po) with po ~ 0.91. Meanwhile, an additional 
pair of imaginary discrete eigenvalues appear for all values of p in [0,7r]. When e — >■ -t-oo, the discrete line soliton 
ipm approaches the slowly- varying function Q with /i = 1, and the interval of real eigenvalues shrinks to {0,po) with 
Po(e) — ^ according to the elliptic 2D NLS equation ([6]) (see and [H, Section 5.9]). 
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FIG. 1: Numerical results for discrete 
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bifurcating from the F point of the dNLS equation ([T]) on a two- 
dimensional lattice. Upper row: profiles of discrete line solitons tprn', middle row: real parts of eigenvalues A of the spectral 
stability problem pO[) versus the transverse wavenumber p; lower row: imaginary parts of eigenvalues A versus p (the shaded 
pink region is the continuous spectrum). Left column: e — 0.1; middle column: e — 1; right column: e — 4. The red dashed 
line in the middle left panel is the leading-order analytical approximation p3p in Theorem [T] 



Next we consider discrete line solitons bifurcating from the X point in the dNLS equation ([T|). At three values 
of e, eigenvalues of the spectral stability problem (1161) for various transverse wavenumbers p in the interval [0, n] are 
presented in Fig. [51 We see that when e — 0.01, a single pair of purely imaginary eigenvalues exist for all values of p 
in (0, tt], in agreement with Theorem [21 These imaginary eigenvalues match the asymptotic formula (jlSp in Theorem 
[21 (lower left panel). When e = 0.2, this pair of imaginary eigenvalues intersect the continuous spectrum (lower middle 
panel). As a consequence, complex eigenvalues appear on the p-interval of 0.97 < p < 1.97 (center panel). When 
e = 4, additional eigenvalues exist. The eigenvalue curves on the left side of the p-interval (right middle and lower 
panels) are the counterparts of similar curves for line solitons in the hyperbolic 2D NLS equation ([7]) (see [10] and 
[29l . Section 5.9]). But the curve of real eigenvalues on the right side of the p-interval (right middle panel) has no 
counterpart in the hyperbolic 2D NLS equation ([7]). These real eigenvalues bifurcate out from the origin inside the 
continuous spectrum. As e — >■ -foo, the eigenvalue curves on the left side of the p-interval shrink toward p = at the 
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asymptotic rate of e Meanwhile, the real-eigenvalue curve on the right side of the p- interval approaches the edge 
point p — n, and its width shrinks at the asymptotic rate of e~^^^. 
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FIG. 2: Numerical results for discrete line solitons (|15|l bifurcating from the X point of tiie dNLS equation ([T]) on a two- 
dimensional lattice. Upper row: profiles of discrete line solitons ^pm', middle row: real parts of eigenvalues A of the spectral 
stability problem (|16l) versus the transverse wavenumber p; lower row: imaginary parts of eigenvalues A versus p (the shaded 
pink region is the continuous spectrum). Left column: e = 0.01; middle column: e = 0.2; right column: e — 4. The red dashed 
line in the lower left panel is the leading-order analytical approximation (|18[) in Theorem [2] 



B. Numerical results for the dNLS equation on a one-dimensional lattice 

Now we consider discrete line solitons ([?T|) in the dNLS equation with k = 1. At three values of e, eigenvalues 
of the spectral stability problem ([22]) for various transverse wavenumbers p are presented in Fig. [3l We see that when 
e = 0.1, a pair of real eigenvalues exist in the interval (0,po), where po ~ 1.53. For p > po, these real eigenvalues 
become purely imaginary. This is in agreement with Theorem |3l Quantitatively, these real and imaginary eigenvalues 
are well approximated by the leading-order asymptotic formula ()24p in Theorem [3l At e = 2, we still have the 
instability band (0,po) with po ~ 1.81. Meanwhile, two additional branches of purely imaginary eigenvalues appear 
over certain p-intervals. When e = 4, the instability band (0,po) has po ~ 1.75, and one additional branch of purely 
imaginary ei genv alues exist over the entire p-axis. When e — > +oo, po{e) — > a/3 according to the elliptic 2D NLS 
equation © [111 . 

Next we consider discrete line solitons (|2T]) in the dNLS equation (IT9l) with k = — 1. At the same values of e, 
eigenvalues of the spectral stability problem (|22p for various transverse wavenumbers p are presented in Fig. 21 We 
see that when e = 0.1, a pair of imaginary eigenvalues intersect the continuous spectrum (lower left panel). As a 
consequence, complex eigenvalues bifurcate out near p — 1/2, in agreement with Theorem 3] (middle left panel). When 
e = 2, additional eigenvalue bifurcations occur (middle column). When e = 4, eigenvalue curves split into two parts. 
The left part is the counterpart of similar curves for line solitons in the hyperbolic 2D NLS equation ([7]) [TO, HI] , while 
the right part is a curve of real eigenvalues at large p. Notice that this eigenvalue structure at e = 4 qualitatively 
resembles that in Fig. [2] (right column) for discrete line solitons bifurcated from the X point in the dNLS equation 
As e — ?> -l-oo, the left part of this structure asymptotically approaches eigenvalue curves for line solitons in the 
hyperbolic 2D NLS equation ([7]). On the other hand, the location of the right real-eigenvalue curve moves to p ^ oo 
at the asymptotic rate of e^^^, and its width shrinks at the asymptotic rate of e~^/^. 

Lastly, we quantitatively compare the numerical complex eigenvalues bifurcating from p — 1/2 with the analytical 
formulae for small e in Theorem 4. For this purpose, we have numerically determined the most unstable complex 
eigenvalue Xmax and its p-location Pmax for each e in the range of < e < 0.3, and the results are displayed in Fig. [S] 
For comparison, the leading-order analytical approximations (|27p for Xmax and Pmax are also plotted in this figure. 
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FIG. 3: Numerical results for discrete line solitons (|2ip in the dNLS equation H19|) on a one-dimensional lattice with «; = 1. 
Upper row: profiles of discrete line solitons i/;™; middle row: real parts of eigenvalues A of the spectral stability problem 



(|22[) versus the transverse wavenumber p; lower row: imaginary parts of eigenvalues A versus p (the shaded pink region is the 
continuous spectrum). Left column: e = 0.1; middle column: e = 2; right column: e = 4. The red dashed lines in the middle 
and lower left panels are the leading-order analytical approximations (|24|l in Theorem [S] 



We can see that the analytical and numerical results closely match each other. 



V. SUMMARY AND DISCUSSION 



In this article, we have analytically determined the transverse stability and instability of line solitons in the dis- 
crete nonlinear Schrodinger equations on one- and two-dimensional lattices in the anti-continuum limit. On a two- 
dimensional lattice, the fundamental line soliton was proved to be transversely stable (unstable) when it bifurcates 
from the X (F) point of the dispersion surface. On a one-dimensional (stripe) lattice, the fundamental line soliton was 
proved to be transversely unstable for both signs of transverse dispersion. In addition to these qualitative stability 
results, we have also derived asymptotic expressions for unstable eigenvalues and compared them with numerical 
results with perfect qualitative and quantitative agreements. 

It is noted that the discrete nonlinear Schrodinger equations are generally used to describe wave dynamics in 
the continuous NLS equations with a deep periodic potential and without inter-band mode coupling. Although the 
analytical results in this article nicely explained many of the numerical results on the transverse stability of line 
solitons in the continuous NLS equations (sol . [3ll | , they cannot explain some other notable facts in the continuous 
models. For instance, our analytical results for the dNLS equation (|19l) on a one-dimensional lattice say that all line 
solitons are transversely unstable, but the numerical results in [31] showed that in the continuous model, line solitons 
near the second Bloch band can be transversely stable. The reason for this discrepancy is that line solitons near the 
second Bloch band contain a strong coupling between the first and second Bloch bands, which is neglected in the 
discrete NLS equation. How to analytically explain the existence of transversely-stable line solitons in the continuous 
NLS equations with a one-dimensional lattice is still an open issue which merits further study. 
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FIG. 4: Numerical results for discrete line solitons (|21[) in the dNLS equation H19p on a one-dimensional lattice with k — — 1. 
Upper row: profiles of discrete line solitons tl^m', middle row: real parts of eigenvalues A of the spectral stability problem 
(|22l) versus the transverse wavenumber p; lower row: imaginary parts of eigenvalues A versus p (the shaded pink region is the 
continuous spectrum). Left column: e = 0.1; middle column: e — 2; right column: e = 4. 
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